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Abstract. We study the behavior of weakly bound clusters and their relation to 
the well-known three-body Efimov states. We adopt a model to describe universal 
behavior of strongly interacting bosonic systems, and we test its validity by reproducing 
predictions of three- and four-body universal states. Then, we extend our study to 
larger systems and identify a series of universal cluster states that can be qualitatively 
interpreted as adding one particle at a time to an Efimov trimcr. The properties of 
these cluster states and their experimental signatures are discussed. 
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The universal behavior of quantum systems has attracted much attention from 
the atomic, condensed-matter, and nuclear physics communities. Theoretical and 
experimental studies strongly suggest that the ultracold behavior of two-component 
Fermi gases is characterized by a single interaction parameter: the scattering length a (cf. 
Ref. [I] and references therein). In contrast, there is no comprehensive understanding of 
universality in bosonic systems with large a. For weakly interacting dilute systems, there 
exists a gaslike state in which the interactions depend on a. However, this state becomes 
unstable against collapse as a increases [2j. At the three-body level (N=3), Efimov 
physics [3] leads to a different concept of universality. In this case, two parameters are 
necessary to describe the low-energy behavior: a and the three-body parameter, k , that 
is related to the energy of an Efimov state at unitarity (\a\ = oo). The experimental 
observation of Efimov phenomena in ultracold gases [I] has reinvigorated the interest in 
universal few-body physics. 

The extension of Efimov physics to larger systems has been an issue of debate. 
Initial studies [5], based on restricting approximations, proposed the existence of an 
"iV-body Efimov effect" for systems with N > 3. This prediction does not agree with 
a full quantum-mechanical treatment of the four-boson problem [6]. Recent four-body 
studies [SI 13 IE] predict that the system is universal, i.e., no new parameter is needed, and 
the low-energy behavior of the system is only described by a and Kq. The experimental 
observation of resonant four-body features [9] , in agreement with theoretical predictions 
of Ref. [6], strengthens our confidence in the understanding of universal four-body 
phenomena. Another series of studies focused on the existence of weakly bound cluster 
states. These studies considered such systems as 4 He and tritium or others that interact 
through realistic van der Waals or other model potentials [TDJ dH [321 EES]- Although 
there are clear differences between these studies, their predictions are qualitatively 
similar, possibly indicating an underlying universal behavior. Furthermore, the studies 
of Ref. [13] suggest that cluster states only depend on two- and three-body physics. 
This opens up the possibility of predicting universal properties of cluster states only in 
terms of a and k [7]. 

Here, we propose a method to explore the properties of bosonic systems with large 
a. We adopt a model Hamiltonian that combines two- and three-body interactions and 
allows the independent control of a and kq. Our model Hamiltonian is based on Efimov 
solutions of the three-body problem and designed to significantly reduce the nonuniversal 
corrections of the lowest energy states. This property makes it particularly suitable for 
studying larger systems for which only the lowest states are computationally accessible. 
Using numerical techniques, such as diffusion Monte Carlo (DMC) and correlated- 
Gaussian (CG) basis-set expansion, we study the properties of clusters. First, we show 
that the model Hamiltonian reproduces the universal behavior of three- and four-body 
systems. Second, we extend our calculations to systems with N < 13 particles, and we 
obtain results consistent with the premise of universality. We conclude that for each 
Efimov state, there is at least one iV-body cluster state. These iV-body cluster states 
are only controlled by two- and three-body physics and inherit many of the properties of 
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the universal three- and four-body states. The cluster states are formed by long-range 
binding and fall into the category of quantum halos [H] since most of their probability 
is in the classically forbidden region. For large a, the energies of different cluster states 
are linearly related, forming Tjon lines [15]. Also, under certain conditions, the iV-body 
cluster is bound, while all smaller clusters are unbound. This behavior manifests the 
Borromean nature of the states [UJ . Our predictions describe states of physical systems 
that fulfill the conditions of universality: large a (|a| larger than any other length scale 
of the interactions) and long-range binding. In experimental studies, these conditions 
might not be fully satisfied by the lowest Efimov family, but can be quantitatively 
fulfilled by excited families whose cluster states are resonances and, therefore, have a 
finite width. 

Model Hamiltonian. — To construct a model Hamiltonian that captures the essence 
of Efimov physics, we consider the analytical solutions of a three-body system with 
zero-range interactions in the hyperspherical framework [16]. In this representation, 
the description of the system reduces to a one-dimensional Schrodinger equation in 
the hyperradius R = \/ (rf 2 + ff 3 + r| 3 ) /3. At unitarity (\a\ = oo), the description of 
Efimov physics is governed by an effective potential U(R) = — (sq + l/4)h 2 /(2mR 2 ), 
where m is the mass of the particles and so ~ 1.00624. The divergence of this effective 
potential at R — > leads to the unphysical Thomas collapse [T7], i.e., an infinite number 
of states with energies approaching — oo. Thomas collapse is an artifact of the zero-range 
interaction being modified when potentials with finite range r are considered. Finite- 
range corrections alter the behavior of U(R) in regions R < rn, producing deviations 
from its universal 1/R 2 behavior and introducing a natural regularization. These finite- 
range effects strongly affect the lowest energy states, introducing important nonuniversal 
corrections. 
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Figure 1. (Color online) Description of the model potential in the hyperspherical 
framework. The vertical line and the shaded region at kqR <C 1 represent the hard- 
wall three-body potential (also included in the inset). The dashed line corresponds to 
the binding energy, and the solid red line to the wave function g^. The dash-dotted 
curve corresponds to the potential U(R). Inset: Comparison of go (solid curve) and go 
(dashed curve) in logarithmic scale. 
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To obtain universal behavior even for the lowest energy states, we adopt a 
regularization scheme that consists of introducing a repulsive three-body interaction 
in the hyperradius. The simplest three-body interaction is a hard wall of size R c , i.e., 
V3 = 00 for R < R c and otherwise. Thus, our hyperradial potential at unitarity is 
U M (R) = -(s 2 + l/A)h 2 /{2mR 2 ) for R > R c , and U M {R) = 00 otherwise. Note that 
U M (R) agrees with U(R) for all R > R c . Thus, the analytical hyperradial solution from 
Ref. |16j . go(nR) = V nRK iso ( y/2nR) , is valid for our model Hamiltonian for all R > R c . 
Here, Kj So is the modified Bessel function of the second kind with imaginary index iso, 
and k is the momentum associated with the energy of the state by E^ b = —K 2 h 2 /m 
(the superscript U refers to unitarity). The hard-wall potential imposes the boundary 
condition go{K n R c ) = that determines the discrete spectrum E^ bn = —K^ff/m. The 
hyperradial wave function of our model Hamiltonian ~% is then g Q (K, n R) for all R > R c 
and otherwise. 

Figure [T] summarizes the basic idea of this regularization scheme. Note that only a 
very small fraction of the wave function go is in the region R < R c (see inset of Fig. [I]), 
suggesting that even the ground state of our model Hamiltonian is, to a large extent, 
universal. For example, the lowest energies follow closely the Efimov-scaling relation 
Eg >n+1 = e- 27T / s0 Eg n . The lowest momenta, k « 0.04622/i^ and k x « 0.0020359/fk, 
imply that the Efimov-scaling relation is obeyed even for the lowest two states within 
0.1%. The long-range binding is another essential signature of universality. The typical 
interparticle distances are much larger than the interaction range R c , which implies that 
the properties of the state are insensitive to nonuniversal short-range physics. These 
properties suggest that even the lowest state can be considered universal and of Efimov 
character. 

To extend the description of universal physics to larger systems, we adopt the 
simplest Hamiltonian that leads to the N-body ground state with long-range binding: 

H = E + E v *(rv) + E (i) 

i i<j i<j<k 

Here is the interparticle distance between particles i and j, and Rijk = 
\J( r ij + r ik + r Jfc)/3 is the three-body hyperradius. The three-body force prevents the 
formation of tightly bound clusters and guarantees long-range binding. This premise is 
confirmed a posteriori when we analyze the pair-correlation function of the cluster states 
and suggests that iV-body forces (with N > 3) are not essential for the description of 
universal physics. 

In Eq. [TJ the ideal combination of interactions would be a zero-range 
pseudopotential, V^(r) = Anah 2 /m5(r)^r, and a hard-wall potential of size R c for 
V3. Numerically, however, these interaction potentials can be hard to deal with; we 
need to replace them with slightly different potentials that are better suited to the 
techniques we apply. For example, the zero-range pseudopotential has to be replaced 
by a finite-range potential. This finite-range potential leads to the same two-body 
physics as the zero-range pseudopotential as long as \a\ ^> r . For CG calculations, 
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we use V^r) = —V^ G e~ r2 ^ 2r o- 1 for DMC calculations, we use a square-well potential: 
Vz(r) = — K) f° r r < r 0) an d otherwise. For DMC calculations, we can use the hard- 
wall three-body potential as the three-body interaction, but for CG calculations we need 
to replace it with V^i?) = V'e~ R2 ^ 2R ^ [6]. For this model, universality is expected in the 
regime \a\ ^> R c ,Tq. To reduce finite-range corrections, we impose the condition that 

r < Re- 
Numerical Calculations. — To study a system with iV < 6, we use a CG basis set 
expansion [TjJJ 02] • In our implementation, the eigenstates of a system are expanded in a 
set of CG basis functions in which the center of mass has been removed, and the relative 
angular momentum is zero. Each basis function is a symmetrized product of Gaussian 
functions, each of which depends on one of the N(N — l)/2 interparticle distances. The 
parameters that characterize the Gaussian functions are selected and optimized using a 
stochastical variational method [19]. This optimization allows us to reach convergence 
of the lowest energy states with a few hundred (thousands) basis functions for N = 3 
(N = 4, 5). 

To extend our calculations to larger systems, we implement a DMC algorithm [20] . 
To within statistical uncertainties, the DMC algorithm provides the exact ground-state 
energy. We introduce importance sampling through a trial wave function ty? to reduce 
the statistical uncertainty. The trial wave function is first optimized using variational 
Monte Carlo (VMC) methods. To describe the weakly bound cluster states, we propose 
the trial wave function: 

* T (ri, r n ) = $(i? T ) J] g(R hhk ) J] /(ry), (2) 

i<j<k i<j 

where / and g are the two- and three-body correlations, $ is a hyperradial correlation, 
Rt is the iV-body hyperradius given by R\ = XX r « ~ R-ca/) 2 , and Rca/ = 
^2,Yi/N is the center-of-mass coordinate. The two-body correlation is the zero- 
energy scattering solution of two particles interacting through V2, i.e., f(r) = (r — 
a) sin(-Ko r )/[ r sin(i^o r o)] for r < r o, and f(r) = 1 — a/r otherwise. Here Kq = mVo/ti 2 . 
The three-body correlation is g(R) = for R < R c , and g(R) = (R — R c )R b exp(—cR 2 ) 
otherwise. Here b and c are variational parameters. The hyperradial correlation is 
= R b ' exp(— c'R 2 ) where b' and d are variational parameters. 

Universal Cluster States. — We extend our studies to larger systems by solving 
the many-body Schrodinger equation of our model Hamiltonian using CG and DMC 
methods. Figure [2] presents the ground-state energies E Nb = E Nb of the bosonic cluster 
state. Considering that these universal clusters are only controlled by two- and three- 
body physics, their energies should obey E^b = n^h 2 /m €n [1/(koo)] = \E^ b \eN[l/(Koa)], 
where the is a universal function that only depends on N. Thus, the numerical 
results as presented in Fig. [2] represent the universal relationship of the bosonic cluster 
energies given by ejy. The symbols represent the Monte Carlo results, while the solid 
curves correspond to other predictions. For N = 3, the curve corresponds to the semi 
analytical prediction of Ref. [16]. For N = 4, the curve corresponds to predictions from 
Ref. [5]. For iV = 5, we carry out CG calculations for the model Hamiltonian. The good 



Weakly Bound Cluster States of Efimov Character 6 

agreement in the N = 3 — 5 comparisons represents strong evidence of the validity of 
our model Hamiltonian. 
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Figure 2. (Color online) Bosonic cluster energies in the negative scattering-length 
region. Symbols correspond to the DMC predictions for R c /r — 4. The numbers at 
the top identify the cluster size, and the vertical dashed lines correspond to l/{na* Nb ). 
For N = 3 and 4, the solid curves correspond to predictions of Ref. [TB] and Ref. [BJ, 
respectively. For N = 5, the solid curve corresponds to CG predictions using our model 
Hamiltonian. Dashed curves are fits to an analytic simple form (see text). 

The predictions from Fig. pimply that there is at least one iV-body state universally 
related to an Efimov state. These states are weakly bound. For example, the trimer 
energies are 10 5 times smaller than the typical two-body interaction energy f?j{mr^). 
Also, their energies follow a smooth and simple dependence on n a that we approximate 
by Em « E^ b (x + cn x bN )/ (1 + cn), where x = (a — a* Nb ) / a, E^ b is the cluster energy at 
unitarity, and a* Nb is the critical scattering length for which the cluster energy reaches 
the zero-energy threshold. E^ b and a* Nb are extracted from the DMC results (see e.g., 
Table [T]) and cn and are fitting parameters. The dashed curves in Fig. [2] shows the 
proposed analytical expression for E eb , E 7b , and E 8b when setting c 6 ~ 0.58, c 7 ~ 0.56, 
eg ~ 0.84 and b& ~ 2.68, bj ~ 2.86, b 8 w 2.50. The simple dependence of E^b suggests 
that, close to unitarity, the different cluster states follow generalized Tjon lines, i.e., their 
energies are linearly correlated to each other and E Nb w + @n Er N -i)b- The Tjon 
lines are verified with the numerics. For example, we obtain E^ f» 2.QQE^ b + 2.QE^. 

The positions of the critical scattering length a* Nb (dashed vertical lines in Fig. [2]) 
reflect the Borromean nature of these states [llj: for a* Nb < a < a1 N _^ b , the A^-body 
cluster is bound but all smaller clusters are unbound. The approximately equal spacing 
of the vertical dashed lines in Fig. [2] suggests that the critical scattering lengths follow a 
simple relationship. Actually, the critical scattering lengths for 3 < iV < 8 can be well 
described by l/(K a* Nb ) ~ 2.3(1) — iV which also suggests that a* Nb — > when N — > oo. 
This prediction implies, in agreement with previous studies [21] . that there exist large 
cluster states for any a < 0. 

Next, we benchmark the energy relations at unitarity. Table [T] presents the DMC 
energies of cluster states up to N = 13 at unitarity. To determine these energies, we 
analyze systems with R c /r = 4, 5, and 8 (which is equivalent to changing k ). The 
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Table 1. Energies at unitarity and scattering-length ratios that characterize weakly 
bound cluster states. The scattering length ratios can be transformed to an absolute 
scale using l/(«o«36) ~ 0.64. 
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0.42(1) 


9 


49.9(6) 
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0.60(1) 
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60.2(6) 
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18.59(5) 


0.71(1) 


11 


70.1(7) 
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27.9(2) 


0.78(1) 
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79.9(3) 
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38.9(3) 


0.82(1) 


13 


88.0(7) 



reported uncertainty is estimated by considering the dependence of the results on R c j r Q 
and the statistical uncertainty. Our predictions agree only qualitatively with those 
of Ref. [33] that considers states more tightly bound and which are more affected by 
nonuniversal corrections. 

The CG studies allow us to analyze the existence of excited bound states. Four- 
body calculations with the model Hamiltonian predict the existence of two states whose 
energies at unitarity are Ef b ~ 4.55E^ b , and E^ b 2 ~ 1.003-E^, in agreement with previous 
studies E]. Next, we extend the study to five particles. First, we determine the 
existence of a five-body state with an energy E^ b rs 10. 4i^, in good agreement with our 
QMC results. Also, we observe an extremely weakly bound excited five-body state. 
This state becomes more strongly bound as nonuniversal corrections become more 
important, i.e., -R c /r — » 0. We conclude that there exists an second state that is bound 
or on the verge of being bound and that, under experimental conditions, nonuniversal 
corrections of the realistic interactions are going to be crucial for ultimately determining 
the existence of the second five-body states. For iV = 6, we can only analyze the ground 
state, but we cannot reach full convergence. However, through an analysis of the energy 
dependence on the basis size and an extrapolation to infinite basis functions, we estimate 
Em ~ 18. 3(3)£/3&, in good agreement with the DMC predictions. 

Resonant enhancement of losses produced by the presence of a bound state at 
the fragmentation energy threshold has been the key signature for identifying universal 
three- and four-body physics jU El |9]. This enhancement can also be a suitable path 
for experimentally observing the universal physics of larger systems. A theoretical 
framework for the analysis of iV-body recombination processes has been recently 
proposed [22]. The critical scattering-length ratios reported in Table [l] complement 
that study by identifying how resonant enhancements produced by different iV-body 
clusters are related to each other. The ratio a* Ab /a\ b is in good agreement with other 
theoretical predictions [6] and experimental results [9]. The ratio a\ h /a\ b predicts that 
for the Cs experiment of Ref. [9] , there should be a five-body resonant feature at about 
0.6 x 410 ao ~ 250ao- The experimental results for the losses show a series of g-wave 
Feshbach resonances in that region that prevent us from identifying a possible five-body 
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Figure 3. (Color online) Pair-correlation distributions for N = 3 — 6. Note that the 
maximum of the pair-correlation increase with TV. For N = 3, 4, circles correspond to 
the DMC predictions, while solid lines correspond to other predictions (see text). For 
N = 5 and 6, the curves correspond to DMC predictions. 



resonant feature. 

Finally, we analyze the structural properties of these cluster states. Figure [3] 
presents the pair-correlation functions, defined as 4irr 2 P(r) = (^f\8(rij — r)\ty), for 
N = 3 — 6 at unitarity. In our DMC implementation, we calculate the structural 
properties, P, using the mixed estimator (P) m ixed = 2(P)dmc ~ {P)vmc [20] . For 
N = 3, we compare the model Hamiltonian predictions calculated with DMC (circles) 
with pair-correlation results extracted from semianalytical results (curve on top of 
circles, indistinguishable from the DMC predictions). For N = 4, we compare the model 
Hamiltonian predictions calculated with DMC (circles) with universal pair-correlation 
results extracted from Ref . [6] . We found excellent agreement in both comparisons. The 
long-range binding of the states can be deduced from the pair-correlation functions which 
show that the mean interparticle distance is much larger than the typical interaction 
distances r and R c . For example, the peaks of the pair correlations are at ~ 7R C (i.e., 
r ~ 30r ) for the parameters used in Fig. [3j 

In conclusion, we have combined a model Hamiltonian with numerical techniques 
to construct a suitable framework for the description of Efimov physics in many-body 
systems. This model description can be extended to study systems with different species 
and spin statistics. Our results present a universal picture for bosons in the a < 
regime and predict a series of cluster states universally related to each other. As a 
consequence of universality, the cluster states inherit the energy and size Efimov-scaling 
relationships that control the universal three- and four-body states. Thus, two clusters 
tied to two consecutive Efimov trimers would have their energies at unitarity related 
by E Nbjn+ i = e~ 2n / s °E Nbjn . These predictions oppose the notion of an "iV-body Efimov 
effect" [5]. Our predictions can be experimentally investigated by analyzing resonant 
losses in ultracold bosonic gases. Alternatively, these cluster states can be created 
by loading bosonic atoms in an optical lattice. The proposal of Ref. [23J of creating 
trimers in optical lattices can be extended to larger systems. In principle, the energies 
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of different cluster states can be measured as functions of a, experimentally reproducing 
the results of Fig. [2] 
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